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Time is a valuable resource and it seems intuitive that longer time should lead to better precision 
in Hamiltonian parameter estimation. However recent studies have put this intuition into question, 
showing longer time may even lead to worse estimation in certain cases. Here we show that the 
intuition can be restored if coherent feedback controls are included. By deriving asymptotically 
optimal feedback controls we present a quantification of the maximal improvement feedback controls 
can provide in Hamiltonian parameter estimation and show a universal time scaling for the precision 
limit of Hamiltonian parameter estimation under the optimal feedback scheme. 


The implementation of quantum technology usually re¬ 
quires a full and precise information about the parame¬ 
ters of system evolution, which makes quantum Hamil¬ 
tonian parameter estimation a crucial problem. An im¬ 
portant task of Hamiltonian parameter estimation is to 
find out the ultimate achievable precision limit with 
given resources and design schemes that attain it. Typ¬ 
ically Hamiltonian parameter estimation is achieved by 
preparing some initial quantum state po and letting it 
evolve under the Hamiltonian H(x), through the evo¬ 
lution p x = U x p 0 Uj., where U x = e ~ lH ( x ) T , the un¬ 
known parameter in the Hamiltonian is imprinted on 
p x , one then can estimate the parameter by measuring 
p x . This problem is well studied in quantum metrol¬ 
ogy when the Hamiltonian is in the multiplication form 
of the parameter H(x) = xH , it is known that in this 
case the optimal strategy is to prepare the initial state 
|A ma „)+|A m i„) ^ where |A max(mm) ) is the eigenvector of 


as 
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H for the maximum(minimum) eigenvalue, the standard 
deviation of the optimal unbiased estimator of x then 
scales as . here n is the number that the process is 

VnJ 


repeated and J = (A ma:r —A m i n ) 2 T 2 is the maximal quan¬ 
tum Fisher information, where A ma x(min) is the rnaxi- 
mum(minimum) eigenvalue of H and T is the time that 
the Hamiltonian acts on initial states)]]. In this case the 
standard deviation of the estimation scales as A which 
shows that longer time always leads to better precision, 
consistent with our intuition. However for general Hamil¬ 
tonian H{x), recent studies have shown different time 
scalings[5j. For example, if a Hamiltonian takes the form 


H(x) = H[cos(a;)(Ji + sin(:r)< 73 ], where o\ 


0 1 ) 
10 ’ 


02 = ^ ^ q* J and <73 = ^ ^ J are Pauli matrices, 

the maximum quantum Fisher information in this case 
is then J = 4sin 2 BT, which oscillates with timej3j[4]. 
Thus longer time may even lead to worse estimation for 
general Hamiltonian, this contradicts with our intuition. 


In this article, we will show that the intuition can be re¬ 
stored when we include feedback controls. By presenting 
an asymptotically optimal feedback scheme for general 
Hamiltonian parameter estimation, we give a quantifica¬ 
tion of the maximal improvement feedback controls can 
provide in Hamiltonian parameter estimation. We show 
that under the optimal scheme the precision limit dis¬ 
plays a universal time scaling ^ which is independent of 
the form of the Hamiltonian. In this article we focus on 
single parameter estimation, generalization to multiple 
parameters is possible but is beyond the scope of this 
article. 

The methods developed previously in 0 SI to com¬ 
pute maximal quantum Fisher information for general 
Hamiltonians are quite invovled and hard to incorporate 
feedback controls. Here we use a tool developed in our 
recent work which is computationally efficient and con¬ 
venient to include feedback controls, which we recapture 
here briefly)!)]. The precision of estimating x from quan¬ 
tum states p x is related to the Bures distance between 
p x and its neighboring states p x +dx\ SHS], 


^BuresiPxi Px+dx) = ~^J(Px)dx , (1) 

here the Bures distance dsures is defined as 

dBures (Pi; P 2 ) = \/ 2 — 2Fb(Pi, p 2 ), (2) 

where Fb(pi, p 2 ) = \J P 1 P 2 P 1 is the fidelity. Thus maxi¬ 
mizing the quantum Fisher information J(p x ) is equiva¬ 
lent to maximizing the square of Bures distance between 
p x and its neighboring states. 

If the evolution is governed by U x = e~ lH ^ x ' T with 
a general Hamiltonian H(x), then p x = U x poUf and 
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Px+dx — U x 4 - dxPO U 3 . _|_ dx ■ thus 

max d 2 Bu res(Px, Px+dx) 

Po 

= uia,xd Bures (U x poU x , U x + dx poU x+dx ) 

Po ^ 

= max[2 - 2F B {U x poUl , U x+dx p 0 Ul +dx )\ ( 3 ) 
= max [2 - 2F B (p 0 , UlU x+dx p 0 Ul +dx U x )) 

= 2- 2mmF B (p 0: U' p 0 U ,t ), 

Po 


where U' = U^.U x+dx . Let e _z6 J be eigenvalues of U ', 
where 9 l - £ (— 7 r, 7 r] for 1 < j < d(here d denotes the di¬ 

mension of U'), which we will call eigen-angles of U', and 
arrange #(^x = 9 B ' > 9 > • • • > 9 d ' = 9^ in in decreas¬ 
ing order. Then min po F B (po, U'poU'^) = cos 


9 U -9 b 


pU 

-t'rv 


if ^max - Cm < We denote C T e{U) = 

for a given unitary operator, then from Eq.Q we get 


G u> — 

8(1 - cos max 0 " mm 


aV 


max J = lim 

Po dx—> 0 

= lim 

dx—t 0 




dx 2 

8[1 - cos C T E{UlU x+dx )\ 


dx 2 


(4) 


The precision limit is then given by 


For Hamiltonians in the multiplication form F[{x) = 
xH, U x = e ~ ixHT , W x U x+dx = e ~ iHTdx , in this case 

FtE B x-{-dx) (^max ^min)-L I dx |, 


Eq.|7| then recovers the well-known formulap] 


Sx > 


y/n( A max ^min )T 


( 8 ) 


For general Hamiltonians, this also provides a straight 
forward way of calculating maximum quantum Fisher in¬ 
formation. We will demonstrate it through an example, 
which will also be used later to show the gain of feedback 
controls. Consider the Hamiltonian H(x) = H[cos(x)cri-l- 
sin(a;)( 73 ], where x is the interested parameter, it can rep¬ 
resent the direction of a magnetic field [31 II] • The Hamil¬ 
tonian can be written compactly as H(x) = B[a(x) • <r], 
where ai(x) = cos(x),ci 2 (x) = 0 , 03 ( 2 ;) = sin(a;). If it 
evolves with time T, then U x = e ~ lH W T = e -iBT[a(x)-o], 
in this case 

U' = UtU x+dx 

_ giBT[d(x)-a] g— iBT[d(x-\-dx)-<r] 


With a simple calculation one can get 

giBT[d(x)-a] g— iBT[a(x-\-dx)-a] ^iB'[a'-a] 


( 10 ) 


Sx > 


n max po J[U x p 0 U£] 


(5) 


1Jtti y/Ul —cos C T e(u1 CC+dx)] rz 

iim dx ^. 0 - 1^1 -V™ 


where n is the number that the measurement procedures 
are repeated. If U x is continuous with respect to x , 
then when dx —> 0, U' = W(x)U(x + dx) —> I, thus 

Cax) ^min 0 and e ° lax ~ 6lmin —► 0. We then get 


max J = lim 

Po dx —>-0 


f) u ' _ o u 

8(1 - cos max min ) 


= lim 

dx—t 0 

= lim 

dx —^0 


16 sin 


dx 2 

n QU 1 _ qU 1 

" max _ min 

4 


dx 2 

(9 U ' - 9 U • ) 2 

\ max min/ 


dx 2 


AC 2 TE (UlU x+dx ) 
dx^O dx 2 


( 6 ) 


which gives the ultimate precision limit 
1 


Sx > 


yjn max Po J 


1 


lim, 


dx —>-0 


2 C T E(UiU*+dx) 

| da: | 


(7) 


here a 1 is a unit vector and 
cos B' = cos 2 (BT) + cos(dx) sin 2 (BT) 

dr 2 

= cos 2 (BT) + (1-—) sin 2 (HT) + 0(dx 3 ) (jq) 

dr 2 

=1 - sin 2 l y BT)% r + 0(dx 3 ). 

Since the eigenvalues of e lB are e ±lB ; we have 

0 mlx = B ' and ^min = ~ B ' ■> thus fl °‘ ax ~ g ° i ° = B'. From 
Eq.Q we then get 

_ cos 

maxJ = lim 8-—- =4sin 2 (HT). 

dx-> o dx 2 

This is consistent with previous studies [310 , however our 
method makes the computation much simpler. 

Next we are going to include feedback controls and 
show that the intuition that time is a valuable resource 
can be restored. We first prove a useful property of 
Cte{U): given two n-dimensional unitary operators Ui 
and U 2 , if Cte{U\) + Cte(U 2 ) < 7r, then Cte(UiU 2 ) < 
Cte(U\) + Cte(U 2 ). To see this, let Ui = e~ lHl and 
U 2 = e~ lH2 where the eigenvalues of Hi and H 2 are in 
[—7T, 7r]. From Thompson’s theorem[21j. there exists two 
unitary operators V\ and V 2 such that 

U\U 2 = e~ iHl e~ iH2 = v 1+v 2 h 2 v2)' (12) 


























3 




FIG. 1: Hamiltonian parameter estimation with feedback con¬ 
trols. 


Denote H 3 = V\H{V 3 + V 2 H 2 V 2 , we have X(H 3 ) -< 
X(H x ) + X(H 2 )&2l, here A (H) denotes the vector whose 
entries are the eigenvalues of H, arranged in decreas¬ 
ing order, i.e., Ai (H) > X 2 (H) > ••• > A n (H) here 
A i(H) denotes an eigenvalue of H at the i-th entry of 
A (H), denotes the majorization relation^], X(H 3 ) -< 
X(Hi) + X(H 2 ) means 

k k k 

E A *(^3)<E A ^^)+E A ^^) (13) 

2= 1 2= 1 2=1 

for 1 < k < n — 1 and 

n n n 

J2k(h 3 ) = J2^ h i) + J2 x ^- 

2=1 2=1 2=1 

From this it is easy to see that 

X 1 (H 3 )<X 1 (H 1 ) + X 2 (H 2 ), (14) 

X n {H 3 )>X 1 (H n ) + X 2 (H n ). (15) 

Thus A!(ff 3 ) - A n (H 3 ) < Ai {H x ) - X n (H x ) + A^) - 
X n {H 2 ). Since Cte( e ~ iH ) = Xll ' H ' , ~ Xn( - H '> ; this shows 

C TE {UiU2) < C TE (U 1 ) + C T e(U 2 ). (16) 

We are now ready to derive optimal feedback controls. 
As shown in Fig [l] under the feedback scheme the evolu¬ 
tion is interspersed with feedback coherent controls. The 
total evolution can be described by 


Umt{x) = U t (x)U x U t (x)U 2 ■ ■ ■ U t (x)U m , 

here U t (x) = with t = U x , U 2 , ■ • ■ , U m are co¬ 

herent controls. Here we assume the controls are evenly 
interspersed, this is mostly for simplicity of notations, 
the analysis can be generalized to uneven case straight¬ 
forwardly. We will first derive optimal controls for the 
case of m = 2 , same strategy works in the general case. 


When m = 2 , U 2 t{x) = Ut(x)UiU t (x)U 2 , then 


U\ t (x)U 2 t{x + dx) 

=U 2 u} {x)u\u\ (x)U t (x + dx)UiU t (x + dx)U 2 
=u\ul{x)U\{Ul{x)U t {x + dx)]U 1 [U t (x)U}(x)\ 
Ut{x + dx)U 2 

=U 2 {U} (x)Ul)[Uj (x)U t (x + dx)](U x U t (x)) 
[Uj(x)U t (x + dx)]U 2 - 


From Eq. © we know finding the maximal quan¬ 
tum Fisher information is equivalent as finding the 
maximum value of CTE[U 2t (x)U 2 t{x + da;)]. Since 
CTE[U 2t (x)U 2t (x + dx)] is determined by the eigenval¬ 
ues of U\ t [x)U 2t (x + dx) and U 2 does not change the 
eigenvalues in this case, so it can be chosen as any uni¬ 
tary. This is reasonable as measurements on the output 
states have already implicitly assumed additional con¬ 
trols. We divide the rest of operators into two parts, 
{U\{x)U\)[U\{x)U t {x + dx)}(lhUt{x)) and U}{x)U t {x + 
dx ), apply the inequality (16) we get 


CTE[U 2 t (x)U 2 t{x + dx)] < C T E{Ul{x)U t {x + dx)] 

+ C TE [(U} {x)U\)[U} (x)U t (x + dx)](UiU t (x))](18) 
= 2 C TE [U}(x)U t (x + dx)], 


where the last equality we used the fact that 
(U}(x)U\)\U}(x)U t (x + dx)](UiU t (x)) has the same 
eigen-angles as U}(x)U t (x + dx). One obvious choice of 
control that saturates the equality is U\ = U}(x), as it 
aligns the eigenvalues of the two parts and the corre¬ 
sponding maximal and minimal eigen-angles add up. In 
this case U 2t (x)U 2t (x + dx) = U 2 [Uj(x)U t (x + dx)] 2 U 2 , 
C T E[U 2t {x)U 2 t{x + dx)] = 2C T e[U} (x)U t (x + da;)]. 

Note that this is not in contradiction with previ¬ 
ous result which showed that feedback controls do not 
help when H(x) = xH]{2\. The reason feedback con¬ 
trols do not help in that case is U{x) commutes with 
U(x+dx), thus (U\x)U\)\U\x)U(x+dx)]{U\U(x)) is al¬ 
ready equal to W(x)U(x + dx) without any controls(i.e., 
when Ui = I), so the eigen-angles are already aligned. 
However for general Hamiltonians, U ( x ) and U(x + dx) 
in general do not commute, in that case feedback controls 
are useful to increase the precision limit. 

This analysis can be extended to general m straight¬ 
forwardly with 

C T E[Ul t (x)U mt (x + dx)] < mC TE [Ul C x)U t (x + dx)], 

(19) 

where the equality can be saturated with the controls 
U\ = U 2 = • • • = U m - 1 = U}{x) and arbitrary U m . In 
practice the true value of x is not known beforehand, thus 
an estimated value x has to be used for the feedback con¬ 
trols, and the controls U x = U 2 = • • ■ = U m - 1 = U] (x) 
need to be updated adaptively. The maximum quantum 
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FIG. 2: Quantum Fisher information for parameter estima¬ 
tion of the Hamiltonian H(x) = f?[cos(i)(Ti + sin(x)<j 3 ] with 
B = 1 ,T = 1 and t = 1/5, i.e., the evolution is inter¬ 
spersed with 5 controls. The controls are taken as Uj ( x ) 
with x = (1 + /3)x and the real value of x is assumed to be 
1. In this case the feedback controls gain as long as /3 is in 
(- 1 . 66 , 1 . 66 ). 

Fisher information is achievable asymptotically when 
x —> x pa na ng ng. In the non-asymptotical regime, 
with a reasonable estimate such strategy can also gain. 
For example with H(x) = B[cos(x)a x + sin(a;)c 73 ] and 
consider the feedback scheme with a total evolution time 
T = mt and let the controls U\ = U 2 = = U](x) 

with x = (l + fi)x, here /3 represents the error of the esti¬ 
mate. In Fig. i> we plotted the Fisher information with 
different /3, it can be seen that such feedback controls 
gain for a broad range of /3. When j3 = 0, the controls 
leads to the maximum quantum Fisher information which 
equals to 4m 2 sin 2 ^, this is asymptotically achievable 
and much higher than the maximum Fisher information 
without feedback controls. When m is sufficiently large, 
sin — = —, the maximum quantum Fisher information 
reaches 4 B 2 T 2 , which means that the ultimate precision 
limit scales as 

Such time scaling is actually universal that holds for 
any H(x) under the optimal feedback scheme, which we 
will now show. Assume the evolution is interspersed with 
to controls, 

Umt{x) = U t {x)U x U t {x)U 2 ■ ■ ■ U t (x)U m , 

here Ut{x) = e,~' lH ^ c ) t with t = Under the op¬ 
timal strategy U x = U 2 = ■ ■ ■ = U m -1 = U}\x ), 
CTE\Um t (x)U m t{x + dx )] attains its maximal value 
mCTE\Ul(x)U t (x + da;)]. If to is taken sufficiently large, 
then t = — can be sufficiently small, in this case 

U}{x)Ut{x + dx) = e iH{x)t e -iH(x+dx)t 

^ [H(x + dx) - H(x)]T (20) 


FIG. 3: Quantum Fisher information for parameter esti¬ 
mation of noisy channel having Kraus operators K\ (x) = 

U(x) and K 2 (x) = ^^ L a 3 U{x) where U(x) = 

exp (-iH(x)t), r] = 0.8 1,/5 , and the Hamiltonian H(x) = 
B[cos(*)cri + sin(ir)cr 3 ] with B = 1, T = 1 and t = 1/5, i.e., 
the evolution is interspersed with 5 controls. The controls are 
taken as U}{x) with x = (1 + /3)x and the real value of x is 
assumed to be 1. In this case the feedback controls gain as 
long as /3 is in (—1.46,1.70). 

thus 

r 1 rrrt/ \tt ( , j m ■ dx) A m ; n (x, dx)]T 

C T E[Ut{x)U t {x + dx)\ = -, 

TO. 

here A max ( m ; n )(a;, dx) denotes the maximum (minimum) 
eigenvalue of H(x + dx) — H(x). Thus when to is suffi¬ 
ciently large 

mC T E[Ul(x)Ut(x + dx)] = [A max (a;, dx) - A min (x,dx)]T, 

which gives the maximal quantum Fisher information 

Tf ^ 3^2 [Amax]^) dx) — A m i n (x,da;)] (o~\\ 

max J(p x ) =4T --, (^l) 

the ultimate precision limit is then given by 

§x > 1 = _ I_ 

VnJ 2 y/n lim dx ^ 0 X ™^ x ’ dx) ~ x X ^ x > dx ) t ’ 

this shows the time scaling for the ultimate precision limit 
is universal. 

Summary and outlook: we derived an asymptotically 
optimal feedback scheme for Hamiltonian parameter es¬ 
timation and showed that under this scheme the ultimate 
precision limit has a universal time scaling, this restored 
the intuition that time is always a valuable resource when 
the evolution is unitary. At the presence of noise, such 
feedback strategy is still useful in increasing the precision 
limit, for example at the presence of dephasing noises as 
shown in Fig]3] the feedback controls still gain for a broad 




























5 


range of /3. Recently it has been shown that for a noisy 
evolution with the Hamiltonian H(x) = xH , quantum 
error correcting techniques can be used to extend the 
time scaling to ^ if the noise is correctable [MHB7] . It 
is straightforward to combine those techniques with the 
optimal feedback controls here to show that universal 
time scaling also holds for arbitrary Hamiltonian under 
correctable noises. The combined method is expected to 
lead to better precision in general, however the optimality 
will generally be lost with a simple concatenation of these 
two techniques(quantum error correcting is not optimal 
in general), finding out the optimal feedback controls at 
the presence of noises will be a future research direction. 
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